Fluctuations of particle ratios and the abundance of hadronic resonances 
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In this letter we will argue that the event-by-event fluctuations of the ratio of positively over 
negatively charged pions provides a measurements of the number of rho and omega mesons right 
after hadronization. This finding can be utilized to put the hypothesis of chemical equilibration in 
relativistic heavy ion collisions to a test. 
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The question to which extent the matter created in rel- 
ativistic heavy ion collisions is equilibrated is central to 
the interpretation of many observables for the existence 
of a new phase of matter. A detailed analysis of the inclu- 
sive single particle yields of several hadronic species has 
led many authors ||-Q| to conjecture that rather early in 
the collision chemical equilibrium has been reached. In- 
deed assuming chemical equilibrium at an early stage of 
the collisions a rather impressive agreement with a large 
body of data can be obtained by adjusting just a few pa- 
rameters, namely the temperature, the baryon chemical 
potential and the strangeness suppression factor (for de- 
tails see e.g. H). However, this analysis has to rely on the 
abundance of final state 'stable' particles and thus has to 
infer the number of most hadronic resonances present in- 
side the system. Some information about the abundance 
of unstable resonances, such as the p and uj meson, may 
be obtained through the observation of the electromag- 
netic decay into dileptons. However, this only provides 
a time integrated yield and thus gives only limited in- 
formation about the abundance of these resonances right 
after hadronization and/or chemical freeze-out. 

In this letter we propose to study the event-by-event 
fluctuations of particle ratios, in particular the ratio 
tt^/tt^ in order to put a strong constraint on the rel- 
ative abundance of some unstable resonance right after 
hadronization/chemical freeze-out. We will show that 
the fluctuations of 7r+/7r^ are quite sensitive to the num- 
ber of the primordial po and lu mesons. In more gen- 
eral terms, the investigation of the event-by-event fluc- 
tuations of particle ratios provides a crucial test of the 
hypothetical chemical equilibration - to see if it also pre- 
dicts two particle correlations correctly in addition to the 
single particle inclusive data. 

The key point to our argument, that the fluctuation 
of the 7r"'"/7r~ is indeed sensitive to the particle num- 
bers at chemical and not at thermal freeze-out, is the 
observation MM that the pion number does not change 
during the course of the evolution of the system through 
the hadronic phase. Typical relaxation times for pion 



number changing processes are of the order of 100 fm/c, 
much longer than the lifetime of the system, which is 
about 10 fm/c. In addition, as we shall argue in more 
detail below, charge exchange processes, which in prin- 
ciple could affect the 7r"''/7r~-fluctuations, lead only to 
small corrections. Finally, considering fluctuations of the 
multiplicity ratio eliminates the effect of volume fluctua- 
tions which are present even with the tightest centrality 
selection. 

Let us now define some notations. We define the fiuc- 
tuation 6Ni by 



N, = (A^,) + 5N, 



(1) 



where (Ni) is the average number of the particle species 
i. Then the variance is given by 

A(iV„ iVO = {6N,SN,) = (iV,2) - {N,f = m, (iV,) (2) 

Here we introduced the notation 

A(iV„ iV,) EE (A, A,) - (A,) (A,) = (SN^j) (3) 

and defined w^ to be the ratio A{Ni,Ni)/ (N^). 

The absence of correlations make the fluctuation in the 
multiplicity very close to the (inclusive) average number 
of particles. Hence, in a classical thermal system, Wi = 1, 
since there are no correlations among the particles. Bose- 
Einstein or Fermi-Dirac statistics introduce correlations 
so that 



w, = 1 ± (nf) I {ui) 



with 

(n) 



cPp 

(27r)3 



n±{p); {n^) 



(fp 

(2^ 



-Xp) 



(4) 

(5) 



where, (+) stands for Bosons, and (— ) stands for 
Fermions H. For the systems of interest here, however, 
the corrections due to quantum statistics are small; for a 
pion gas at a temperature of 170 MeV, w^ = 1.13 |^. 

In general however, wt will not be equal to 1 ± 
(nf) / {ui) due to additional correlations introduced by 
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interactions and resonances. This will be discussed in 
detail below. 

Given the above notation, the fluctuation of the ratio 
is given by p2[ 






' A(JVi,iVi) ^ A(iV2,7V2) 

^ A(JVi,iV2) 
(A^i) {N2) 



(6) 



The last term in Eq. (^ takes into account correlations 
between the particles of type 1 and type 2. This term 
will be important if both particle types originate from 
the decay of one and the same resonance. For example, 
in case of the 7r+/7r~ratio, the po, uj etc. contribute to 
these correlations. Also, volume fluctuations contribute 
here. 

(i) Volume fluctuations: Even though data are of- 
ten selected according to some centrality trigger, the im- 
pact parameter and thus the volume of the created sys- 
tem, still fluctuates considerably. Assuming, that the 
particle abundance scales linearly with the volume, vol- 
ume fluctuations translate directly into fluctuations of 
the particle number. 

However, by considering ratios of particles these fluc- 
tuations cancel to leading order. This can be seen as 
follows. Note that we can rewrite Eq. (pi) as 



A{Ri2,Rl2) 

{Ruf 



5Ni 6N2 



,(^1) (^2), 



(7) 



Assuming that the volume fluctuation separates from the 
density fluctuation, we write 



N = n,^,i{V)+SV)il + 6n/n,,,) 



where n„ 



SN 



(N) I (1/). Then to the first order. 



6V_ 

(V) 



6V_ 

{V) 



{N) 



(8) 



(9) 



where SnN = (V) dn is the number fluctuation due to 
the density fluctuation. Clearly, the volume fluctuation 
part will be cancelled in Eq. (Q) and hence A{Ni, Nj) in 
Eq. (pf) can be simply replaced with 



A^{N„Nj) = {SnN,6,,Nj 



(10) 



From now on, unless otherwise signified, we will omit the 
subscript n from An{Ni,Nj). 

Let us now turn to the discussion of the density fluc- 
tuations. In the physical system we consider, the den- 
sity fluctuation is mainly due to the thermal fluctua- 
tion. As already mentioned above, in absence of any res- 
onances/interactions the thermal fluctuations A(iVi, A^i) 
are simply given by Eq. (0), with wi slightly differ- 
ent from unity due to quantum statistics. Furthermore 



in a thermal system the correlation term vanishes, i.e. 
A{N,N,) = 5^j A{N,,Nj) since (iVi N2) = (A^i) {N2) in 
that case. This changes, however, once interactions, in 
particular resonances, are present in the system. 

(ii) Effect of resonances: A fundamental assump- 
tion of the statistical model is that at the chemical freeze- 
out time, all the particles including resonances are in 
thermal and chemical equilibrium. The expansion af- 
terwards breaks the equilibrium. However, as discussed 
above the total number of n^ and tt" given by 



{N,)^{N,)^+Y,{R)T{^^)l 



(11) 



remains constant from this time on. Here the subscript 
T on {N.i)T and {R)t denotes the average number of par- 
ticles and resonances at the freeze-out time and {ni)j^ is 
the average number of the particle type i produced by 
the decay of a single resonance R. 

The presence of resonances which decay into the par- 
ticles of interest affects the fluctuations of each individ- 
ual particle {A{Ni,Ni)) M. Resonances decaying into 
both particle species of interest also affect the correla- 
tion term {A{N.i,Nj)). A single resonance contribution 
to A{Ni,Nj) is given by ||ll| 



Ar{N,,N,) 

(R) 



= {n^nj)jj + {WR - 1) {ni)j^ {nj)j^ (12) 



where we defined 



{niU-i) 



3/R 



J2 b?infUnf\ 



rebranches 



and 



{■n-i. 



rebranches 



b?inf)r 



(13) 



(14) 



Here the index r runs over all branches, b^ is the branch- 
ing ratio of the r-th branch, and {nf)r represents the 
number of i-particles produced in that decay mode. At 
the chemical freeze-out, the system is in equilibrium and 
hence the number of different particle species are uncor- 
related. In the final state where all the the resonances 
have decayed, the correlation is given by 

A{N,,N,)=wf{N^^6,j 

+ Y. ^^)t {{n^n3)R + (wl ~ 1) (n,)^ (n,)^) (15) 

R 

where wf {Ni)^ denotes the part of the variance due only 
to the statistical fluctuations at the chemical freeze-out 
time. 

Putting everything together we get for the fluctuations 
of the ratio 



A{Ri2,Rl2) 

{Ruf 



1 



{N2) 



(D 



11 



D22 - 2D 



12; 



(16) 



with 



Z?22 = F2 

D12 = ^ {nin2)j^ 



{R)t 



(17) 
(18) 
(19) 



where 



F,= 







= il + {wf~ l)ri + J2 {{^')r - ("»)fl) 






(20) 



and we defined r.i = {Ni)rp/{Ni). Here we regarded 
(w|J — 1) to be neghgibly smaU. For a typical resonance 
of mn ~ 1 GeV, (w^ — 1) < 10~^ assuming the tem- 
perature of 170 MeV. (In the numerical results presented 
below, these terms, though small, are included.). Note 
that we have factored out 1/ {N2) to separate out the 
explicit dependence on the system size. 

Before we turn to the practical applications, a few com- 
ments are in order at this point. First, the effect of the 
correlations introduced by the resonances should be most 
visible when (A^i) ~ (A'2) since the branching fraction nf- 
should enter with about the same weight. In this case, 
a resonance decaying always into a pair of particles "1" 
and "2" contributes about equally to Dn, D22 and D12 
and hence contributes neghgibly to (-D11+-D22 — 2_Di2). 
On the other hand, the presence of such resonances does 
influence the total number of particle "2", (7V2). Hence, 
for instance, p° and to will always reduce the fluctuation 
in 7r+/7r~ compared to the statistical fluctuation. Sec- 
ond, when (A^2) ^ (-^1)1 as in the K to tt ratio, the 
fluctuation is dominated by the less abundant particle 
type and the resonances feeding into it. 

After this general formulation of the problem, let us 
now turn to the calculation of the negative to positive 
pion ratio. The formalism developed above can be easily 
applied to the case of 7r+/7r~fluctuations by setting 



Ni 



No 



Typically, ^7r"'")/^7r~) ~ 1. In table (|) we show the 
most important contributions from hadronic resonances. 
This calculation includes mesons and baryons up to 
771 ~ 1.5 GeV as listed in the particle data book. Weakly 
decaying strange particles are regarded as stable, but let- 
ting them decay changes the main result very little. The 
values of temperature, baryon chemical potential and the 
strangeness chemical potential are the same as those in 
H, i.e., T = 170 MeV, fib = 270 MeV and ^, = 74 MeV. 
As already pointed out, the p'^ and uj contribute 
about 50 % of the correlations. Furthermore, the 



correlation-term 2Z)_| is seventy percent of the in- 
dividual contributions -D++ or D This is a siz- 
able correction which should be visible in experi- 
ment. Furthermore, since only very few resonances 
have decay channels with more than one charged pion 
of the the same charge, (77±)„ — {n±)j^ ~ to a 
good approximation. Also, {w'^ — l)r± ~ 4 % at 
T — 170 MeV. Hence, the fluctuations of the indi- 
vidual pion contributions are very close to the statis- 



tical limit of D 



++ 



!?__ 



1. Thus the resonances 



contribute predominantly to the correlation term D^ 

TABLE I. Contributions from diflerent hadrons to the 
fluctuations of the tt^/tt" ratio. Contributions are given 
in fraction of the total result. 



Particle 



Others 



0.32 



0.00 



0.02 



0.08 



0.08 



0.00 



0.07 



0.44 



D_ 



0.00 



0.32 



0.02 



0.00 



0.08 



0.08 



0.07 



0.44 






0.00 



0.00 



0.06 



0.00 



0.24 



0.00 



0.21 



0.55 



0.31 



0.00 



0.02 



0.09 



0.09 



0.00 



0.08 



0.43 



0.00 



0.31 



0.02 



0.00 



0.09 



0.09 



0.08 



0.44 



TABLE II. Total values 



D\ 



Dt 



2DX 



Values 1.09 



1.09 



0.76 



Using (7r+) = (ti"^) — (tt^ -|- tt^) /2 and (7r+ + tt^) = 
220 from the recent NA49 results ||l^ on event-by-event 
fluctuations of the transverse momentum and the K/n- 
ratio, we obtain for the total fluctuation tt^/tt^ ratio 



A(i?+_,i?+_ 



{R+-r 



- = al_= 0.0128 



(21) 



so that a^ = 0.113 A more useful quantity, to consider, 

however, is the ratio of the above total fluctuations over 
the purely statistical value. The latter can be possibly 
obtained in experiment by the analysis of mixed events. 
Possible artiflcial contributions due to experimental un- 
certainties, such as particle identification etc. should can- 
cel to a large extent in this ratio, and thus a comparison 
with theory becomes more meaningful. The value for the 
statistical fluctuation for us is simply given by 



1 



1 



(^+) {n-) 



(22) 



since in mixed events all correlations, even the quantum 
statistical ones, are absent, i.e. w+ = ?«_ = 1. Our result 
for this ratio is 



{afl'Y 



0.70 



(23) 



Note that this ratio is independent of (tt" ) . Thus the cor- 
relations introduced by the presence of resonances reduce 
the fluctuations by 30 %. Looking at table (Q), the most 
important contribution of the correlation come from the 
p and Lo meson. Or in other words, the measurement 
of the fluctuations of the tt^/tt^ provides a strong con- 
straint on the initial number of p'^ and w-mesons. 

Doing the same analysis for the K /-k ratio, where pre- 
liminary data exist po[ , we find 



'K/t^ 



(„stat \2 



1.04 



(24) 



in good agreement with the data. Note, however that 
our value for the fluctuation itself, (7^/^ — 0.17 differs 
from the experimental value of cr|^f = 0.23 indicating 
the effect of additional fluctuations from particle identi- 
fications etc. ]1^. The reason (see |11 for details), is that 
in case of the K /i: ratio, the fluctuation is dominated by 
the contribution from the kaon, which is largest due to 
the small abundance. 

Let us close by discussing some possible caveats. First, 
there is the question of acceptance cuts. Clearly, one 
should not consider the particle ratio of the full An accep- 
tance. In this case charge conservation will impose severe 
constraints and reduce the fluctuations of the 7r+/7r~ ra- 
tio. However, as long as a limited acceptance in, say, ra- 
pidity is considered, the constraint from charge conserva- 
tion is minimal and our assumption of a grand-canonical 
ensemble are well justified. Limited acceptance on the 
other hand may reduce the effect of resonances on the cor- 
relation term, as some of the decay products may end up 
outside the acceptance. In order to estimate this effect we 
have performed a Monte-Carlo study. Given the rapidity 
distribution of charged particles for a Pb+Pb collision at 
SPS-energies ||l^, a rapidity window of Ay = 1 changes 
the above results by less then one percent. On the other 
hand this window covers only a fraction of the observed 
rapidity distribution, so that constraints from charge con- 
servation are negligible. Finally there are the charge ex- 
change reactions such as 7r+-|-7r^ ^^ Tr^-t-Tr". These re- 
actions in principle could change the 7r+/7r^-ratio in a 
given event. However, detailed balance requires the net 
change to be close to zero. In addition, for a given event, 
these reactions infiuence not the difference between 7r+ 
and TT^ but only their sum. Thus, again we expect only 
small corrections to the above result since the sum is 
large compared to the expected changes. 

The reaction with baryons, notably tt~+p <-> n^+n 
might be more effective. But again, detailed balance and 
the fact that the nucleon/pion ratios is so small should 



make these corrections insignificant. A detailed quanti- 
tative investigation of these corrections will be presented 
in0. 

In conclusion, we propose to study the event-by-event 
fiuctuations of the 7r+/7r~-ratio in relativistic heavy ion 
collisions. This measurement will provide important in- 
formation about the abundance of short lived resonances 
right after hadronization and/or chemical freeze-out. It 
will further impose a strong test on the validity of the 
chemical equilibration hypothesis in these reactions. If 
chemical equilibrium is reached with the values for tem- 
perature and chemical potential extracted from the single 
particle distributions, we predict that the fluctuations of 
the 7r+/7r~ -ratio should be about 70 % of the statistical 
ones. It would be also of great interest to study these 
fluctuations in proton-proton and peripheral heavy ion 
collisions, where the particle abundances also seem to in- 
dicate chemical equilibrium. 

In addition, any value of the a^ /ct^*°' significantly 

larger than 1 cannot be explained with a simple hadronic 
gas picture and thus would indicate new physics. One 
possible scenario might be the Quark Gluon Plasma bub- 
ble formation |1J]. Thus the 7r+/7r~ ratio as a function 
of Et may serve as an alternative signal for the QGP. 
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